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Abstract. 

This paper provides a general theory to characterize and construct a decoherence- 
free (DF) dynamics for infinite dimensional (i.e., continuous- variable) linear open 
quantum systems. The main idea is that, based on the Heisenberg picture of 
the dynamics rather than the commonly-taken Schrodinger picture, the notions of 
controllability and observability in control theory are employed to characterize a DF 
dynamics. A particularly useful result is a general if and only if condition for a linear 
system to have a DF component; this condition is used to demonstrate how to actually 
construct a DF dynamics in some specific systems. It is also shown that, as in the finite 
dimensional (i.e., discrete- variable) case, we are able to do coherent manipulation and 
preservation of a state of a DF dynamics. 
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1. Introduction 

The theory of decoherence-free subspaces (DF subspaces or DFSs) [H E] provides a very 
promising basis for realizing various key technologies in quantum information sciences 
such as quantum computation [21 SI |5l El [Tj [HI |9], memory [lOl E], communication 
[121 [131 [H] , ^-iid methodology [I5]. It also should be mentioned that some experimental 
demonstrations have been reported [HI [TTl [TH [19]. The usefulness of the theory of 
DPS lies in the fact that, although in reality any quantum system is an open system 
interacting with surrounding environment, the existence of a DF component in a given 
open system implies the realizability of an ideal quantum information processing based 
on the corresponding closed system. 

Let us now review the basic idea of DFSs. Our interest is in a finite-dimensional 
open quantum system interacting with an environment composed of finite number 
of boson field channels. Especially when the interaction is instantaneous and thus a 
Markovian approximation can be taken, the system's density matrix p evolves in time 
with the following master equation: 

rl" ™ 1 1 

= -i[Hrp\ + E {UpLl - -Llhp - -pLlU), (1) 

k=l 

where H is the system Hamiltonian and Lj. represents the coupling of the system to 
the fc-th environment channel. Due to the second term, the time-evolution of p is not 
unitary, implying that in general coherent manipulation or preservation of an arbitrary 
state of an open system is impossible. Nevertheless, if the system has a specific structure, 
which is usually related to a symmetric structure of the dynamics, the system contains 
a component that is not affected by the environment. More specifically, if the system 
matrices H and Lk satisfy specific algebraic conditions [H [2], then in a certain basis 
the system state is represented as p = diagjpDF, 0} (block diagonal matrix) with /5df 
obeying a unitary evolution, while the decohered (D) component p = diag{0,pD} is 
still subjected to a non-unitary dynamics; that is, the dynamics of these states are 
respectively given by 

dppF 

dt 
dpD _ 

dt 

k=l 

where Huf is the Hamiltonian that governs the unitary time-evolution of pdf- Therefore, 
devising H-^y appropriately, we can carry out some desirable unitary gate operations 
even in an open system. 

Based on the above discussion, we here come up with the following question; is 
it possible to construct a DF dynamics of the form (|2l) for an infinite-dimensional 
system? This question is equivalent to asking the applicability of the above theory 
for discrete variable (DV) systems to continuous variable (CV) systems [201 [21]. Recent 
rapid progress of quantum information sciences with infinite-dimensional/CV systems 



-i[ifDF,PDF], (2) 

-z[^',Pd] + {l'JbL',* - -L[*L',p^ - -PbL',*L',), (3) 



3 



strongly indicates the need of a general theory of DFSs for those systems. However, the 
decomposition p = diag{pDF; 0} makes sense only when p is a finite dimensional matrix; 
in contrast, it is very unclear how to make such a partitioning of an infinite-dimensional 
operator. A particular situation where the theory for finite dimensional systems is still 
applicable to a given infinite dimensional system can happen when we represent that 
system in terms of a countable basis such as the number basis {|^)}n=o,i,2,... [21 l22l [23] . 
Actually then the cutting-line between the DF and the D spaces, if it exists, becomes 
clear. Note that this approach is rigorously valid only when the state of the DPS can be 
completely represented by a finite set of basis vectors; but even a very simple coherent 
state is not contained in this class of states. Therefore, for such a case, we need to 
truncate the set of basis vectors and obtain an approximated DF state. 

This paper provides a basic framework for examining a general linear quantum 
system having a DF component, without any approximation. A linear system is an 
infinite dimensional open system that describes the dynamics of, for instance, optical 
modes of a light field [SUES], position of an opto-mechanical oscillator [SUET], vibration 
mode of a trapped particle [28], collective spin component of a large atomic ensemble 
[291 [301 [31] , oi' more generally canonical conjugate pairs of a harmonic oscillator network 
[32l [33l [Ml [35] , which can all serve as quantum-information devices. The key ideas are 
twofold. As mentioned above, in the infinite dimensional case, there is no general 
way of partitioning a density operator into the DF and the D parts explicitly in the 
corresponding infinite-dimensional Hilbert space. The first main idea is that, therefore, 
we focus on the set of physical observables that completely characterizes the system of 
interest, rather than its state, and divide them into two classes corresponding to DF and 
D parts. For this purpose, we represent the system in the Heisenberg picture; in this 
case the system dynamics is not described by the master equation ([1]), but by a quantum 
stochastic differential equation (QSDE) [211 [251 [361 [371 [38] . the QSDE formalism, the 
dynamical equation explicitly contains the input noise brought from the environment 
and also the output equation is obtained in an explicit form. These features are critical, 
because we are then able to naturally define a DF observable as follows; i.e., it is an 
observable such that the sub-dynamics composed of DF observables is not affected by the 
input noise and that the behavior of those observables cannot be seen from the output. 
The second main idea is that we employ the notions of controllability and observability in 
control theory [391 HO] , in order to characterize such an isolated sub-dynamics of a given 
system. (Note that the term "observability" is conventionally used in control theory, 
and it does not mean a physical quantity.) In particular, in the linear case, it is well 
known in control theory that such controllability and observability properties can be 
straightforwardly captured by simple matrix algebras that only require us to compute 
the rank of certain matrices. As a result, based on the QSDE representation of a given 
linear open quantum system, we obtain a simple necessary and sufficient condition for 
that system to have a DF component; this result will be given in Section 3. 

Once a DF dynamics is constructed, then we should be interested in the stability 
of the dynamics. This is a property that autonomously erases the correlation between 
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the DF and the D parts, and it has been extensively investigated for finite dimensional 
systems [HI |l2l US]. In our case, by limiting the state to a Gaussian state 011 H5] . 
a simple if and only if criterion that guarantees the stability of a DF linear system is 
obtained, which will be discussed in Section 4. 

The remaining part of the paper is devoted to the study of some concrete examples 
of a DF linear system, which is divided into two topics as follows. 

(i) (Section 5) The first topic deals with coherent manipulation and preservation 
of a state of a DF linear system. More specifically, it is demonstrated that, in a specific 
two-mode linear system having a one-mode DF component, arbitrary linear unitary gate 
operation and state preservation are possible under relatively mild assumptions. Also 
it will be shown, for the same system, that preserving a state of the DF dynamics is 
equivalent to preserving a certain entangled state among the whole two-mode system. 
In addition, a stability property is examined, showing its importance in realizing robust 
state manipulation and preservation with the DF system. 

(ii) (Section 6) The second is about how to engineer a DF dynamics for a given open 
linear system. We actually find in two specific examples that, by appropriately devising 
an auxiliary system coupled to the original one, it is possible to embed a physically 
meaningful dynamics in the DF system. Therefore, in addition to the practical merit 
that a DF system can be used for coherent quantum information processing, a DF 
system is useful in the sense that it can be used for simulating an ideal closed quantum 
system even in a realistic open environment. 

We use the following notations: for a matrix A = (aij), the symbols A"^ , A'^ , and 
A^ represent its Hermitian conjugate, transpose, and elementwise complex conjugate of 
A, i.e., A"!" = (a*j), = (aji), and A^ = {a*j) = (^4"'')^, respectively. For a matrix of 
operators, A = (aij), we use the same notation, in which case d*j denotes the adjoint 
to CLij. In denotes the n x n identity matrix. 3? and denote the real and imaginary 
parts, respectively. O is a zero matrix with appropriate dimension. 

2. Preliminaries 

2.1. The QSDE 

Let us consider an open system interacting with some environment channels, which 
are assumed to be independent vacuum fields for simplicity. The time-evolution of an 
observable of this system can be described by a QSDE. Let aj(t) be the annihilation 
operator of the i-th vacuum field and take the Markovian approximation; i.e., we 
assume that ai{t) instantaneously interacts with the system and satisfies the canonical 
commutation relation (OCR) [aj(s), a*(t)] = 5ij5{t — s). Then, the field annihilation 
process Ai{t) = ai{s)ds satisfies the following quantum Ito rule [211 [251 ESI [371 [38] : 

dAidA* = 5ijdt, dAidAj = dA*dA* = dA*dAj = 0. 

The system-field coupling in the time interval [t, t + dt) is described by the unitary 
operator U{t + dt,t) = exp[^^{LidA* — L*dAi)], where Lj is a system operator 
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representing the coupling with the i-th vacuum field. A system observable X{0) evolves 
in the Heisenberg picture to X{t) = U{t)*X{0)U{t), where U{t) is the unitary operator 
from time to t, which is constructed from the relation U{t + dt) = U{t + dt,t)U{t). 
Then the QSDE of X{t) is given by 

dxit) = {t[Hit),xit)] + J2 {L:{t)x{t)U{t) - -L*{t)U{t)x{t) - -^t)L:{t)U{t)))dt 

i=l 

m 

+ J2 {[X{t),Ut)]dA* - [X{t),L*{t)]dA.). (4) 

i=l 

Here, we have added a system Hamiltonian H and defined H(t) = U(t)*HU(t) and 
Li{t) = tj{t)*LiU{t). The master equation ([T]) of the unconditional density operator of 
the system, is obtained through the time evolution of the mean value of X{t) and 
the relation {X{t)) = Tr(X(0)/3(t)). The change of the field operator is also obtained 
explicitly; the output field := U{t)*Ai(t)U(t) after the interaction satisfies 

dA^'^^t) = Li{t)dt + dAi{t). (5) 
2.2. Quantum linear systems 

In this paper, we consider a general multi-mode boson system with canonical conjugate 
pairs {qi,Pi). Assume here that the system is composed of n subsystems, and then define 
the vector of operators x := {qi,Pi, ■ ■ ■ , qn^VnV ■ The CCR [qi^Pj] = iSij leads to 

xx^ - {xx^y = iT.n, S„ := diag{S, . . . , S}, S := ^ ^ ^ . (6) 

(S„ is a 2n X 2n block diagonal matrix.) We are interested in a system whose Hamiltonian 
and couphng operator are respectively given by 

H = x'^Gx/2, Lk = cjx, 
where G = G M2nx2n ^ Q2n |25j_ Also let us define the quadrature pairs of 

the noise process: 

4 = (A + i*)/v^, h = {A.,-A*)/V2i, (7) 

and collect them into a single vector as VV = (Qi, A, • • • , Qm-, PmY ■ Then, from the 
QSDE (jlj), the vector of system variables x{t) = {qi(t) , pi(t) , . . . , qn(t) , pn(t))~^ , where 
qi{t) = U{tyqiU(t) and pi{t) = U(t)*piU(t), satisfies the following linear equation: 

dx{t) = Ax{t)dt + S„C^S„c/W(t). (8) 

The coefficient matrices are given by A = + C'^S„C'/2) G R2"x2" and 

C := ^/2(g?(ci), 53(ci), . . . , 3fJ(cJ, %cjy G M^™^^". 

Moreover, corresponding to the quadrature representation of the input field ([7]), let us 
define the output process as 

gout _ ^^out ^ A°'''*)/V2, A°"* = (^r - ir*)/v^^, 
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and collect them as VV°"* = (Qr*> • • -^QT^Pm^V- Then, from equation (P we 
have 

rfyi;°"t(t) = Cx{t)dt + dyv{t). (9) 

Note that all the 2m elements of the vector >V°"* cannot be measured simultaneously; 
the output equation corresponding to a linear measurement (e.g., Homodyne detection) 
is of the form 

dY{t) := MdVV°"*(t) = MCx{t)dt + MdWit), 

where M is a m x 2m real matrix satisfying MS^M^ = O, guaranteeing that all the 
elements of Y{t) commute with each other as well as with those of Y{s) at different 
time s. To denote the system variable, in what follows we often omit the time index t 
and simply use x\ to avoid confusion, the initial value is explicitly written as x(0). 

Now let (x) be the mean vector and V = {AxAx'^ + {AxAx~'')~^) /2, Ax = x—{x) be 
the covariance matrix, where the mean is taken elementwise. Note that the uncertainty 
relation V + i^njl > holds. For the linear system ([8]), the time-evolutions of (x) and 
V are respectively given by 

d{x)/dt = A{x), dV/dt = AV + VA^ + D, (10) 

where D = STC'TC'S„/2. 

2.3. Gaussian systems 

A quantum Gaussian state can be characterized by only the mean vector (x) and the 
covariance matrix V [HJ HS]. More specifically, its Wigner function is identical to the 
Gaussian probabilistic distribution with mean (x) and covariance V. A notable property 
of a linear system is that it preserves the Gaussianity of the state; that is, if the initial 
state of the system ([8]) is Gaussian, then, for all later time t the state is Gaussian with 
mean {x{t)) and covariance V{t), whose time evolutions are given by equation ( fTOll . A 
steady state of this Gaussian system exists only when A is a Hurwitz matrix, i.e., all the 
eigenvalues of A have negative real parts. If it exists, the mean vector is {x{oo)) = 0, 
and the covariance matrix ^(oo) is given by the unique solution to the algebraic matrix 
equation AV (oo) + V (oo) A^ + D = O. Note that a non-Gaussian state can be generated 
in a linear system. 

2.4. Controllability and observability 

Let us here consider the following general linear system, which is not necessarily limited 
to a quantum system: 

X = Ax + Bu, y = Cx, (11) 

where x is the system variable of dimension n, and u and y are an input and an output 
with certain dimensions, respectively. Note that in this case u represents any input such 
as a disturbing noise and a tunable control signal. Now define the controllability matrix 
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C = {B, AB, . . . , A^~^B) and assume that rank(C) = n'. (Note n' < n in generaL) 
Then, there exists a hnear coordinate transformation x x' = {x[~^ ,x'2^)~^ such that 
the system ( fTTj) is transformed to 

^21 ^22 ) ) \ ^2 j 

where x'^ is of dimension n — n' . This equation shows that x'^ is not affected by m, 
hence it is called uncontrollable with respect to (w.r.t.) u. Of course, if n' = n, or 
equivalently if C is of full rank, then all the elements of x is affected by u and in 
this case the system is called controllable. Next, let us define the observability matrix 
O = (C^, A'^C^, (A^)"-!^^)^ and assume that rank(C) = n". Then, there exists 
a linear coordinate transformation that transforms equation ( ITTj) to a linear system of 
the following form: 




A" A" 

^11 ^12 

O A'i, 




y = (o, c'^ 




where x" is of dimension n — n". Clearly, x" cannot be seen from the output y, hence 
x'l is called unobservable w.r.t. y. Note that, as remarked in Section 1, "observable" is 
a term conventionally used in control theory; we should not interpret that x'( is not a 
physical quantity. 

The notions of controllability and observability play the fundamental roles in 
characterizing various important properties of a system, such as stabilizability of the 
system via the input or capability of constructing an estimator that continuously 
monitors the system variables. See e.g. [39l HQ] for full description of those theories. 



3. Decoherence-free mode of a linear quantum system 

This section provides our first main result; that is, based on the Heisenberg picture 
description of a given linear quantum system, a general if and only if condition for that 
system to have a DF component will be presented. 

As briefly discussed in Section 1, in the finite-dimensional DFS can be 

well defined as follows [H [2]. An open quantum system having a DFS is a system 
whose density matrix p, which is subjected to the master equation ([T]), is represented 
in a specific basis as p = diagjpDF, 0} with p^F governed by the unitary dynamics 
(E]). The DFS is a Hilbert space spanned by the basis of pdf- Within the QSDE 
formalism, this means that any observable fully characterized by these basis vectors is 
completely isolated from the environment, thus it is not affected by the input VV(t) and 
the output field ([5]) does not contain any information about that observable. That is, 
a variable (observable) that is decoherence free should be defined as a variable that is 
uncontrollable w.r.t. VV(t) and unobservable w.r.t. >V°^*(t) in the sense described in 
Section 12.41 In this paper, we focus on the infinite dimensional linear system given by 
equations (IHl) and (IHl); a formal definition is given below. 
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Definition 3.1 Suppose that the vector of observables x subjected to the linear quantum 
system ^ and ^ can be decomposed, in a specific basis, to x = {xj)p,x^)~^ where xdf 
is uncontrollable w.r.t. W and unobservable w.r.t. VV°"* while xd is complement to xbf- 
Then, if xdp satisfies the CCR it is called the decoherence-free (DF) mode. 

The purpose here is to fully characterize the linear system having a DF mode. As 
in the classical case, the following controllability matrix C and the observability matrix 
O will play a key role in deriving such a characterization; 

O = {C^, A^C^, . . . , (A^)^"-!^^)^. 
The following lemma will be useful: 

Lemma 3.1 v e Ker(O) <^ e Ker(C^) 

Proof: First suppose v G Ker((9). This means that CA^v = 0, VA; > 0, which is 
further equivalent to C(S„G')^f = 0, Vfc > 0, due to A = + (5^2^(5/2). This 

condition yields 

for all A; > 0, hence S„f G Ker(C^). The reverse direction is readily obtained. ■ 
The following theorem is our first main result: 

Theorem 3.1 The linear system ([^j and ^ has a DF mode if and only if 

Ker{0) n Ker(OS„) ^ 0. (12) 

In this case, there always exists a matrix Ti satisfying range(Ti) = Ker(O) nKer(OE„) 
and T^J^nTi = S^, and the DF mode is given by xdf '■= T^x. 

Proof: We start to show that equation ( IT2|) is an iff condition for the existence of an 
uncontrollable and unobservable mode. First, the only if part is proven as follows. Now 
Ker(0) 7^ is assumed, hence we take the maximum set of vectors spanning Ker((9) 
and denote them by Vi, {i = !,...,£). Then, from the assumption, a subset of them, 
Vi {i = 1, . . . ,i; i < i), is contained in Range(C)'^, where the subscript c denotes the 
complement of the set. Note here that 

Range(C) = {Range(C) n Ker(C)} U {Range(C) n Range(C^)}. 

Clearly, Vi is orthogonal to all the elements of Range(C) fl Range(C^). Also, Vi can 
be always taken so that it is orthogonal to all the elements of Range(C) fl Ker((9). 
As a result, Vi is contained in Ker(C^). Now, from Lemma |3.H {i = !,...,£) 

span Ker(C'''), hence Vi {i = !,...,£) can be represented by a linear combination of 
them: Vi = Ylk=i ^ik^nVk-, which leads to S^Wj = — Yl^k=i ^ik'^k- Then, by the definition 
Ovi = 0, we have OT^nVi = 0; these two equations imply equation f|T2l) . 

The if part is straightforward to prove; from the assumption there exists a vector 
V satisfying v G Ker((!?) and S„f G Ker((9), but from Lemma [3.11 the latter condition 
means v G Ker(C^). Thus, the system has an uncontrollable and unobservable mode. 
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Now we prove, by a constructive method, that there always exists a DF mode under 
the condition ( !T2|) . First, let us take a vector vi G Ker((9) fl Ker((9S„), which readily 
implies that E„t;i is also contained in this space. Note that vi and S„f i are orthogonal. 
Next we take a vector V2 G Ker(C) nKer(CE„) that is orthogonal to both vi and 
Then, as before, T,nV2 is also an element of this space, and further, it is orthogonal to vi, 
S„t>i, and V2- Repeating the same procedure, we construct Ti = (t>i, S„t>i, . . . ,vi, S„t>£); 
this matrix satisfies T^Ti = l2e and S„Ti = TiS^. Hence, we have Tj^S^Ti = S^, 
implying that xdf = 7"]^ a; satisfies the CCR ([6]). ■ 

Below we give an explicit form of the dynamics of the DF mode, by actually 
constructing the linear coordinate transformation from a; to x' = {x]^p,x]^)~^ . First, 
a vector ui is taken in such a way that it is orthogonal to all the column vectors of 
Ti; then is also orthogonal to all the column vectors of Ti and ui. Second, we 

take a vector U2 that is orthogonal to both ui and T^nUi in addition to the condition 
mJTi = 0; repeating this procedure, we construct T2 = . . . , Un-i, S„'u„_£) that 

satisfies T/S„T2 = O and T2''S„T2 = S„_£. By construction, is complement to 
the uncontrollable and unobservable mode, hence we can set xd Definition 
13.11 From the proof of Theorem 13.11 together with the above discussion, the matrix 
T = (Ti,T2) yields an explicit form of the dynamics of x' = T~^x. Noting that T is 
orthogonal and symplectic, which leads to T^S^ = S^Tj^ and T2^S„ = T^n-iTj , we find 
that equation (IHl) is transformed to 




= IV Am, T2) dt+ IV SnC's^rfw 




S^Ti' GTi O 



dW, (13) 



where the relation CTi = O is used. Also, T2 ATi = O follows from the condition 
CA^Ti = 0, Wk > 0, which implies there exists a matrix P satisfying ATi = TiP. 
Moreover, the output equation ([9]) is represented as 




^y^out _ ^^^^^ ^1.. + _ CT2Xj,dt + dW. (14) 



From equations (fT3|l and ([HD, xdf is apparently a DF mode that is completely isolated 
from the environmental noise fields. In particular, the equation dx-op/dt = S^T/GTiXdf 
means that the DF mode obeys the unitary time evolution governed by the Hamiltonian 

^DF = xdf(0)^G'df£df(0)/2, Gdf = T^GT,. (15) 

That is, the unitary gate operation U-DF(t) = exp(— it^DF) can be carried out coherently 
with the DF mode xdf- 

Lastly in this section, let us consider a practical question; given a system- 
environment interaction, what kind of system Hamiltonian is allowed in order for the 
system to have a DF mode? More precisely, for a given C, if the condition f|T2l) is 
satisfied when G = O and a suitable matrix Ti has been assigned, can we add a nonzero 
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matrix G such that equation f ll2l) still holds and the same DF mode xdf = T^'x exists? 
Here a simple characterization of such G is given: 

Proposition 3.1 Suppose that, for a given C, there exists a matrix Ti satisfying 
Range(Ti) = Ker(C') fl Ker(C'S„) and S„Ti = TiS^. Then, z/Range(Ti) is invariant 
under a linear transformation G, the system specified by this G and the same C has a 
DF mode T^x. 

Proof: The condition means that there exists a matrix Q satisfying GTi = TiT^jQ, 
which is equivalent to GTi = T,^TiQ due to the assumption on Ti; thus (S„G')Ti = TiQ 
holds. Then, since CTi = O, we have C{IlnG)^Ti = O for all k > 0. Again 
from the relation T^nTi = TiS^, this implies C'(S„G')'^S„Ti = O, VA; > as well. 
Consequently, it follows from a similar calculation shown in the proof of Lemma 13.11 
that OTi = OUnTi = is satisfied; this means T^x is a DF mode of the system. ■ 

This proposition further leads to a convenient result in a particular case: 

Corollary 3.1 In addition to the assumptions given in Proposition \3.1\ suppose 
Ker(C') = Ker(C'S„). Then, the system specified by G and the same C has a DF 
mode T^x, if and only if CGTi = O. 

Proof: Because of Range(Ti) = Ker(C'), the condition C'(E„G)'=Ti = O, VA; > 
holds if and only if (E„G')Ti = TiQ, 3Q. Then, due to the same reason, we equivalently 
obtain CGTi =0. ■ 



4. Stability of the decoherence-free mode 

In this section, we assume that the system's initial state is Gaussian, hence the state 
is always Gaussian and is characterized by only the mean vector and the covariance 
matrix, as described in Section 12. 3[ Based on this assumption, we here study a stability 
property of the dynamics of the DF mode xdf- Note that in the finite-dimensional 
case some useful results guaranteeing similar stability properties of a DFS have been 
obtained in 011 ill |13] . 

In order to coherently manipulate a quantum state of £df, i-e. Pdf, it has to be 
separated from that of xd; actually, if xdf is correlated with x d, this means that pdf can 
be affected from the environment. This separability of the states of DF and D modes 
is, in Gaussian case, equivalent to the fact that the covariance matrix V of the whole 
system has to be expressed a.s V = diagjl^F, Vd}, where Vdf and Vy) correspond to the 
covariance matrices of xbf and xd, respectively. Now, from equation ( ITOl) . the covariance 
matrix of the system (fT3|) changes in time with the following matrix differential equation: 

dV f Ai O \ f Aj O \ ( O 0\ 

where Ai = S^Cdf, ^2 = T^ATa, and D2 = T^J:'^C^C^nT2/2. Hence, if the 
solution of equation f[T^ takes a block diagonal form V{t) = diagjVDFl^), ^(^)} in 
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a long time limit without respect to the initial value V{0), this means that the DF 
dynamics is robust in the sense that an unwanted correlation between the DF and the 
D modes autonomously decreases and finally vanishes. The following theorem provides 
a convenient criterion for the system to have this desirable property {eig{A) denotes any 
eigenvalue of a matrix A). 

Theorem 4.1 3ft[eig(S^G'DF) + eig(A2)] < if and only if in a long time limit the 
solution of equation [Td\) takes a form V{t) = diag{l^F(i), Vb(t)}, where Vdf(^) and 
Volt) correspond to the covariance matrices ofx-£,p{t) andxY){t), respectively. Moreover, 
if Gbf > 0, then the above iff condition is simplified to 3f?[eig(A2)] < 0. 

Proof: When partitioning the matrix variable V a.s V = (Vi, V2; Vg^, V3), the 
dynamics of V2 is given by 

dV2/dt = A1V2 + V^Aj. 

We take a vector form of this matrix differential equation. For this purpose, let V2 be 
a collection of the row vectors of V2, which is a real 4:i{n — £)-dimensional row vector. 
Then t>2 obeys the dynamics of the form dv2/dt = {Ai (g) l2(n-i) + he ® A2)v2. Now note 
that any eigenvalue of A = Ai^l2[n-e) + he^A2 is equal to eig(Ai) + eig(A2). Therefore, 
we have f2 — ?■ as t — > 00, if and only if the real parts of sum of any eigenvalues of 
Ai = S^Gdf and A2 are strictly negative. 

To show the second part of the theorem, let us consider the eigen-equation 
^eGoFd = ^9 with g the eigenvector and A the corresponding eigenvalue. From this we 
have (7^(GdfS^G'df)5' = ^g'^G-opg, which immediately leads to (A + A*)(7^G'df5' = 0. If 
Gbf > 0, this can be further expressed as (A + A*)||a/Gdf5'P = 0, implying A + A* = 
or GBPg = 0, thus 3fJ[eig(E£G'DF)] = 0. ■ 

We remark that a specific assumption G^f > is significant from the viewpoint of 
state manipulation; recall that the unitary gate operation with the DF mode is given 
by UDFit) = exp(— zt^DF), which is defined below equation ( |T5i) . Then, it was shown 
in [35] that, if Gdf > 0, for any e > and Xi < there always exists T2 > such 
that ||t/DF(Ti) — Udf{t2)\\ < e with || • || an appropriate operator norm. This result 
means that a unitary operation requiring propagation with negative time can be always 
realized by a physically realizable unitary operation with positive time; this condition 
is indeed necessary to implement a desirable unitary gate operation via a sequence of 
unitaries generated from a set of available Hamiltonians. 

5. Coherent state manipulation and preservation with DF mode 

Coherent state manipulation and preservation are certainly the primary applications of 
a DFS, as demonstrated extensively in the finite dimensional case. This section presents 
two simple examples to show how a DF mode can be used in order to achieve the same 
goals with an infinite dimensional linear system. 
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Figure 1. Two trapped particles in a two-sided optical cavity. 



5.1. Particles trapped via dissipative coupling 

Let us consider a system composed of two particles such as ions trapped in a two-sided 
optical cavity depicted in figure [1] The coupling operator between the particles and 
the outer environment field is L = y/K,{ai + 02), where ai = {qi + ipi)/\/2; qi and pi 
represent the position and momentum operators of the i-th particle. This dissipation 
effect described by L stems from the coupling of the particles to the single-mode cavity 
field that immediately decays and can be adiabatically eliminated; the parameter k 
represents the decay rate of this dissipation. The system dynamics ([8]) and the output 
equation ([9]) are now given by 

dx = Axdt + BdW, rf>V°"* = Cxdt + dW, 

We begin with the assumption that the system does not have its own Hamiltonian 
(i.e., G = 0) and then determine a DF mode. Now C = 1/ k/2{1, i, 1, i), thus when 
G = we have 

o = c = v-J' ° ' °], OE. = cE. = y;;f ' ' ° 

^\^oioiy' ^ y-1 -1 J 

Therefore, Vi = (1, 0, -1, 0)^/^2 and S2t;i = (0, 1, 0, -1)^/^2 span the kernels of both 
matrices; i.e., span{fi, S2fi} = Ker((9) = Ker((9S2). Hence, from Theorem 13.11 the 
system has a DF mode. In particular, the procedure shown in the proof of Theorem 13. II 
yields Ti = (t>i, ^2^1), thus the DF mode is given by 

xdf = X = — = . . . 

V2 \P1~P2 J 

Here let us discuss adding a nonzero matrix G that does not change the above- 
obtained DF mode. Now Ker((9) = Ker((9E2), hence from Corollary 13. G is required 
to satisfy CGTi = O. Then it is immediate to see that G has to be of the form 
G = (Gi, G2; G2, Gi) where Gi and G2 are both real 2x2 symmetric matrices. Therefore, 
the DF mode is driven by the Hamiltonian ( TTSl) with Gdf = T^GTi = Gi — G2. This 
means that, if both of the particles are governed by the Hamiltonian represented by 
Gi and this Hamiltonian can be arbitrarily chosen, then we are allowed to do arbitrary 
linear gate operation on the state of xdf- In particular, by setting Gdf = O, we can 
preserve any state of the DF mode. 
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Next, to evaluate the stability of the DF mode, let us explicitly construct the 
dynamics of both xdf and xd, according to the procedure shown above equation f|T3|) . 
To obtain the matrix T2, we chose Ui = (1, 0, 1, 0)~^ /\/2 as a basis vector orthogonal to 
Vi and T.2V1. Then T,2Ui = (0, 1, 0, 1)^/ \/2 is orthogonal to Vi, E2W1, and Ui, implying 
T2 = (ui, S2U1). Hence, the linear coordinate transformation matrix is obtained as 




T=iT,,T2), 

The transformed variable x' = T~^x then satisfies 
dx' = A'x'dt + B'dW, = C'xdt + dW, 

£(6-1 - G2) O 



(17) 



A' 



O 



E{Gi + G2) - 



B' 




C' 



As expected, the first CCR pair of x', i.e., £df, is neither affected by the incoming noise 
field nor observed from the output field. Now assume that the state is Gaussian. Then, 
by taking large /t so that the condition of Proposition 14. II is satisfied, the stability of the 
DF mode is guaranteed. That is, the off diagonal block matrix of the covariance matrix 
V' = {Ax'Ax'~^ + {Ax'Ax'^)~^)/2 converges to zero as t — )■ 00, thus in the long time 
limit the states of DF and D modes are always separated. Hence, with this Gaussian 
system, arbitrary linear gate operation and state preservation can be carried out in a 
robust manner, under relatively mild assumptions. 

Lastly, based on the assumption that the state is Gaussian, let us examine a state 
preserved in the whole system, under the condition Gdf = Gi — G2 = O. In fact, with 
this condition the DF mode xdf does not change in time, and hence the whole state 
of X is preserved if the D mode is at steady state. This condition is achieved by again 
taking large k so that the coefficient matrix of the drift term of Xb is Hurwitz; in this 
case, the D mode has a steady Gaussian state with covariance matrix 12/2. As a result, 
the covariance matrix of the whole state of x' is V'{0) = diag{VDF(0), -^2/2}, where here 
the time is reset to t = 0. The covariance matrix of x is then given by 



V{0) = TV'{(})T^ 



I2/2 + Vdf{0) I2/2-VMO) 
/2/2-Vdf(0) /2/2 + Vdf(0) 

If the initial state of the DF mode can be set to a squeezed state with covariance matrix 
Vi)f{0) = diag{e72,e-72}, then ^(0) takes the form 



(18) 



and, as mentioned above, this does not change in time. Thus the state preserved is a pure 
entangled Gaussian state; actually, the purity for this state is P = 1/ ^yl6det{V{0)) = 1, 
and the logarithmic negativity |16], which is a convenient measure for entanglement, is 
Ej\f = r/2 > 0. This fact implies that, in general, preservation of a nontrivial state of 





( l + e"- 




1 


1 + e-^ 


1 - er '' 


4 




1 + e'' 






1 + e-^ / 
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a DF mode can lead to that of an entangled state of the whole system, if the stability 
of the DF mode is guaranteed. In other words, a system having a stable DF mode can 
be utilized for preserving an entangled state among the whole system, in addition to 
preservation of a state of the DF mode. 

Remark: It was found in a different and specific framework [47j that, for the same 
system, the Gaussian state with covariance matrix f|T8|) is preserved . 



5.2. Particles trapped via dispersive coupling 



The second example shown here deals with the two particles trapped in an optical cavity 
as before, but with different coupling L = yJ~Kj2{c[\ + ^2), where again qi is the position 
operator of the z-th particle. This coupling realizes a continuous-time measurement of 
^1 + ^2; see e.g. [IHl |49] for how to actually construct this kind of position monitoring 
scheme. Now C = a/ fi;/2(l, 0, 1, 0), hence, when G = 0, we have 




OS2 




and thus 



Ker((9) = span 



( 1 





-1 



Ker(C»S2) = spanj 



V / 

/ 1 \ 
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/o\ 
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The intersection of these two spaces is the same as that obtained in the previous 
example, thus we have the same coordinate transformation matrix ( !T7|) . But of course 
the dynamical equation of x' = T~^x differs from the previous one: 

dx = A'xdt + B'dW, dW"""^ = C'x'dt + dW, 

^ ^ B' = -V2-J ^ 
O S(Gi + G2) \E22 ^ 

where En = diag{l,0} and E22 = diag{0, 1}. Note that = (^i + Q2)/V^ is not 
affected by the noisy environment while it appears in the output equation. This means 
that we can extract information about q'2 without disturbing it; that is, ^2 is a QND 
observable. 

Let us now discuss the state manipulation and preservation with this system. As in 
the previous example, if two particles have the same and arbitrary Hamiltonian specified 
by Gi, this means that any linear gate operation on the DF mode is possible. However, 
this system does not satisfy the condition of Theorem 14. 11 implying that the system does 
not have capability to decouple the DF and the D modes once they become correlated. 
In this sense, the state manipulation with this DF mode is not robust. Next, to see if 
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Figure 2. Schematics of (a) an opto- mechanical osciUator coupled to an optical cavity 
field and (b) two opto-mechanical oscillators coupled to a common optical cavity field. 

the system can preserve a state of the whole system, let us again set the initial state 
to the Gaussian state with covariance matrix ( !T8|) . A notable difference to the previous 
case is that now the state is a time-varying one with covariance matrix 

V^(O) ^ Vit) = ViO) + ] ( ^'l^' ] . 

The purity of this Gaussian state is P = 1/ ^ylQdet{V{t)) = 1/ v^l + Kt and thus 
converges to zero as t — )■ oo. That is, any useful state does not remain alive in the 
long time limit. As a result, due to the lack of the stability property, this system offers 
neither desirable state manipulation nor preservation. 

6. Engineering a DF linear system 

This section deals with the topic of how to engineer a DF dynamics, when a certain 
linear open system is given to us. More precisely, the problem is not how to find a DF 
component in that system, but rather how to synthesize an auxiliary system so that the 
extended system has a desirable DF mode. 

The main purpose of this problem is to devise a DF system for coherent state 
manipulation or preservation as discussed in the previous section, but we here provide 
another important perspective regarding this subject. Needless to say, an ideal closed 
system representation is very strong in describing and exploring the nature of the system 
of interest; however, in reality, any quantum system is essentially an open system 
interacting with surrounding environment. Therefore, generating a subsystem that is 
exactly isolated from environment is of vital importance. 

6.1. Opto-mechanical oscillator with DF mode 

The first example is an opto-mechanical oscillator depicted in figure [2] (a); the system 
is composed of a two-sided optical cavity where one of the mirrors is movable and acts 
as an oscillator. Let qi and pi be the oscillator's position and momentum operators, 
respectively. Also define q2 = {0,2 + 0,2) /V^ and p2 = (0-2 — a2)/y/2i with 0,2 the 
annihilation operator of the cavity mode. The oscillator is driven by the Hamiltonian 
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H = mu'^qf+pl/m where m and u represent the mass and the resonant frequency of the 
oscillator, and further, it is coupled to the cavity mode through radiation pressure. The 
cavity mode couples to an outer coherent light field at the other mirror in a standard 
dissipative manner. As a result the dynamical equation of the composite system is given 
by 

dx = Axdt + bfdt + BdW, = Cxdt + dW, 



A 








1/m 








\ 




—mu'^ 
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and h = (0, 1,0,0)^. Here, 7 denotes the coupling strength between the oscillator and 
the cavity, and k the decay rate of the cavity field. Also f(t) is a real- valued classical 
force acting on the oscillator, which is sometimes taken as an unknown force to be 
estimated; but here we assume f(t) = 0. See e.g. ^26j for a more detailed description of 
the system. 

It is immediate to see that O is of full rank, and thus the system does not have 
a DF mode. Hence, as posed in the beginning of this section, we try to engineer an 
appropriate extended system so that it has a DF mode. In particular, we follow the 
idea of Tsang [50]; an auxiliary system with single mode (^3,^3) is prepared and directly 
coupled to the cavity mode with strength g, leading that the extended system of variable 



A, 



is given by 

: AeXedt + BedW, 



dW 



out 



CeXedt + dW, 
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where the auxiliary system is assumed to have the Hamiltonian [—i'q^ + np'^)/2. Here O2 
denotes the 2x2 zero matrix. The problem is to find a set of parameters {g, fi, v) such 
that the above extended system has a DF modeQ. First, we need that the observability 
matrix O composed of A^. and C*g is not of full rank; this requirement yields /iz/ = — w^. 
Next, the iff condition given in Theorem 13. II imposes us to have ^ = 1/m. As a result. 



the parameters must satisfy fi = 1/m and u 



-muj 



implying that the Hamiltonian 



of the auxiliary system is exactly the same as that of the original oscillator. Although 
we can fabricate this auxiliary system using only some optical devices such as a DPA, 
a natural situation is that another opto-mechanical oscillator with the same mass and 

I The problem considered in [50j is to find the same parameter set (g, /i, v) so that the input noise field 
Q{t) does not appear in the output field P°^^{t), for the purpose of enhancing the sensing performance 
of the incoming unknown force f{t). 
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Figure 3. Three particles trapped in a ring-type cavity. 



resonant frequency serves as the auxiliary system. In this case, the (6, 3) entry of 
should be taken as hg, which yet does not change the above selection of the parameters. 
The resulting extended system is shown in figure [2] (b). Note that the coupling strength 
g can be chosen arbitrarily; this strikingly differs from ^50j where g is also required to 
satisfy a certain condition. 

The dynamical equations of the DF and the D modes are explicitly obtained as 
follows. With the parameters determined above, the procedure shown in Section 3 
yields the linear coordinate transformation matrix T = (Ti,T2) as 



Ti 
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Y 







\ 
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ih 


O2 


\ 




O2 




T - ^ 




O2 


ih 
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gh 
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where 7' = i + g^- Hence, the DF mode is given by 



xdf 



1 

Y 



gqi - 1% 
gpi - 7P3 



Together with the D mode, the transformed dynamics of x'^ ■ 

C'xjt + dW, 



T~^Xe is given by 



dx'^ = 


A'Xdt + B'JW, 






/ 1/m 






K = 


—muj'^ 
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BL 



Oo 
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-b: 



where now the parameter 7 in A is replaced by 7' = a/Y~+^- ^ notable point is that 
the DF mode is governed by the same Hamiltonian as that of the original oscillator. 
This means that an ideal closed opto-mechanical oscillator can be in fact implemented 
in a realistic open system. 



6. 2. Trapped particles with DF mode 

Let us again consider the system of two particles discussed in Section 15. 1[ which is now 
subjected to the following Hamiltonian: 

^ = IY1 ^^p'j + ^^'^D + ^^(^1 - ^2)', (19) 

where m and u are the mass and the resonant frequency of both of the particles. The 
particles interact with each other through a harmonic potential with strength k, as 
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indicated in the second term of equation 0191] . The corresponding G matrix is of the 
form G = (Gi, 6*2; ^2, Gi), thus from the result obtained in Section |5TT] this system 
has a DF mode. However, this single-mode DF system cannot, of course, simulate the 
behavior of the closed two-mode system driven by the Hamiltonian ( IT9|) . Therefore, let 
us try to engineer a suitable extended system having a two-mode DF dynamics with 
which the closed dynamics driven by the Hamiltonian (|T9|) is realized. In other words, 
we want to perform a coherent simulation of the closed system with Hamiltonian ( !T9|) . 
in an open quantum system. 

Here we consider an extended system composed of three particles trapped in a 
ring- type cavity depicted in figure El The Hamiltonian is assumed to be 

He = H + \i-^pl + rnuj'^(il) + \H(12 - %f + \Hq^ - qi)\ 

where kj denotes the coupling strength between the particles and uj' the resonant 
frequency of the auxiliary particle; these parameters will be appropriately chosen 
later. The particles couple to the common single-mode cavity field that can be 
adiabatically eliminated, which as a result leads to the system-bath interaction described 
by Le = y/K,{ai + 02 + 03). The system matrices in this case are given by 

f + k + ks ~k -fcg \ 

—k u"^ + k + k2 —k2 , 

-ks -k2 w'^ + k2 + ks J 

and C = Here, for simple notation, the matrices are represented in the 

basis of {qi,q2,q3,pi,p2,P3)~^ this notation is used throughout this subsection. Also, 
now m = 1 is assumed without loss of generality. 

If G = O, the linear coordinate transformation matrix T = (Ti, 

yields the DF mode xdf = x and the D mode xd = T2 x. Here the problem is to find 
a condition of G 7^ O so that this xdf is indeed a DF mode even when the Hamiltonian 
specified by G is taken into account. Now Corollary 13.11 is applicable, hence G must 
satisfy CGTi = 0, and this yields co' = uj. In this case, it can be actually verified that 
xdf is the DF mode driven by the Hamiltonian H-^y = a^DpG^Dpa^DF/S with 




G — I ^ ^ I , Gq 




UJ 



+ 2k + {k2 + A;3)/2 V3(/C3 - k2)/2 



Gdp - T,GT, - I ^^^^ _ ^^^^2 + 3{k2 + h)/2 1 ® 

We are particularly interested in the question whether the coupled particles with 
Hamiltonian ( !T9l) can be simulated with the isolated DF mode xdf = x. Then, 
because the two particles are identical, the condition 2k + (^2 + kz)/2 = 3{k2 + k3)/2 
is necessary. Further, choosing the coupling strength as k2 = \/3fc and k^ = (2 — \/3)k, 
we end up with the expression 



+ V^k + (3 - V3)k -(3 - V3)A; 



Gdf-I ' Vq^ ' '/Q^7/ "^''^ ,,2 ^'-''^ 1®^2. 



-(3 - V^)k ^2 ^ ^ (^3 _ ^3)^ 
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The corresponding Hamiltonian driving the DF mode is thus given by 

i=i,2 

where now xdf = ('?i5 ^2)P'i)P2)^- Recall that m = 1 is assumed. We can then prove 
that Gdf > and also that A2 = T2AT2 is Hurwitz. Therefore, from Theorem 14. H the 
DF mode with the above Hamiltonian is stable. That is, although the parameters differ 
a bit from those of the target system, we can realize the isolated two-mode coupled 
particles having the same structured Hamiltonian as that of the original open linear 
system, with the guarantee of stability. 

7. Conclusion 

In this paper, we have developed a basic theory of general DF linear quantum systems, 
shown some applications, and demonstrated with some specific examples how to actually 
construct a DF dynamics. The theory contains a general iff condition for a linear system 
to have a DF mode, which fully utilizes the notions of controllability and observability 
in control theory. Although, as expected from the literature in the finite-dimensional 
case, there should be a number of applications of DF linear systems, in this paper we 
have studied only very basic coherent manipulation and preservation of a state of a 
single-mode DF dynamics. Rather what was focused on is the idea of synthesis of an 
open linear system having a desirable DF component, to maintain the usefulness of the 
general results obtained in the former part of the paper. Several applications of the DF 
linear systems theory will be presented soon. 
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